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Abstract 

The degenerate Blume-Emery-Griffiths (DBEG) model has recently been introduced 
in the study of martensitic transformation problems. This model has the same 
Hamiltonian as the standard Blume-Emery-Griffiths (BEG) model but, to take into 
account vibrational effects on the martensitic transition, it is assumed that the states 
S = have a degeneracy p (p = 1 corresponds to the usual BEG model). In some 
materials, the transition would be better described by a disordered DBEG model; 
further, the inclusion of disorder in the DBEG model may be relevant in the study 
of shape memory alloys. From the theoretical point of view, it would be interesting 
to study the consequence of conflicting effects: the parameter p, which increases 
the first-order phase-transition region, and disorder in the crystal field, which tends 
to diminish this region in three dimensions. In order to study this competition in 
high-dimensional systems, we apply a mean-field approximation: it is then possible 
to determine the critical behavior of the random DBEG model for any value of the 
interaction parameters. Finally, we comment on (preliminary) results obtained for 
a two-dimensional system, where the randomness in the crystal-field has a more 
drastic effect, when compared to the three-dimensional model. 
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1 Introduction 



The hamiltonian of the degenerate Blume-Emery-Griffiths (DBEG) model 
reads: H = -J 'E<o> SiSj - KJ2<i,j>SfS] + AJ2iSf, where J is the ex- 
change constant, K is the biquadratic interaction, A is a crystal field, the 
first two sums are over nearest-neighbor pairs on a lattice, the third sum is 
over all sites, Si = ±1,0, and the 5 = states have a degeneracy p [1]. This 
parameter mimics the effects of vibrational degrees of freedom on marten- 
sitic transitions. The DBEG model was studied, for a particular value of the 
crystal-field, within a mean-field approximation and Monte Carlo simulation: 
the effect of p is to shrink the ferromagnetic phase and to increase the region 
where the transition is of first-order [1]. 

Nevertheless, some materials were found to be better described by a disordered 
DBEG model, which may be relevant in the description of conventional shape 
memory alloys as well [1]. From the theoretical point of view, on the other 
hand, it has been shown that randomness may have drastic consequences 
on multicritical behavior. In two dimensions, for instance, any infinitesimal 
amount of disorder supresses non-symmetry-breaking first-order phase tran- 
sitions and replaces symmetry-breaking first-order phase transitions by con- 
tinuous ones. The effect of disorder on high- dimensional systems is different: 
first-order phase transitions disappear only at a finite amount of randomness 
[2]. This behavior has been observed in a large variety of models, either with 
random "field" or with random "interaction" [2-4]. 

Thus, in the disordered DBEG model two ingredients have opposite effects: 
while the degeneracy factor p increases the region where the transition is of 
first-order, disorder tends to supress this region. This competition may have 
interesting effects on the phase diagrams, apart from potential experimental 
relevance. In order to study this competition on high-dimensional systems, we 
apply a mean-field approximation to the DBEG model in a random crystal- 
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field, which follows the probability distribution V(Ai) = r S(Ai + A) + (1 — 
r) S(Ai — A). It is worthy stressing that, if the interaction parameters, J 
and K, were chosen to be random, instead of the crystal field, the overall 
consequences on the phase diagram would be the same [2] . 

2 Formalism and Results 

Our approximation is based on the Gibbs inequality for the free energy [5]: 
F < TrpTi + (l//3)Trpln p, where F is the exact free energy, Ji is the DBEG 
hamiltonian and p is an exactly solvable density matrix. We chose p as the 
most general single-site density matrix for the ferromagnetic (positive J and 
K) system. Since this procedure is usual, we will not discuss it in detail. 

Let us first comment on some general results we obtain, for any K/J. The 
transition at zero temperature is not affected by the parameter p and takes 
place at (A/z J) = (K/J + 1)(1 + r)/2; for smaller values of A/z J the mag- 
netization equals 1, while m = r for A/zJ > (K/J + 1)(1 + r)/2. 

For A/z J = oo, the random field DBEG model, with the field distribution 
given by the equation in the previous page, is equivalent to the quenched site- 
diluted spin-1/2 Ising model. The states S — ±1 represent magnetic sites and 
the states 5 = represent non-magnetic impurities. Therefore, only for r >r c 
(where r c is the lattice-dependent value of r where an infinite cluster of present 
sites is formed for the first time, when r is increased from zero) the critical 
line between 02 and D (see figures) should extend to A/ J = oo. For r < r c , 
the "tail" of continuous transition would touch the zero temperature axis at 
a finite value of A/z J. The reason we find this "tail" extending to A/ J = oo 
for any r ^ is that the present approximation is equivalent to a model with 
infinite-range interactions (see Ref. [3]). 

The assymptotic value of kT/zJ for A/z J = oo equals r. This result comes 
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from the mapping cited in the last paragraph: the degeneracy of the S = 
states and the biquadratic interaction K play no role in the dynamics of the 
S = ±1 states. If the distribution V(Ai) = r<J(A i ) + (l-r)<J(A i -A) is used, a 
different result will hold. The limit A/zJ = oo is now equivalent to the spin-1 
site-diluted Ising model; therefore, K and p are relevant and the assymptotic 
value of kT/zJ depends on both parameters. 

Since for high values of K/J the critical behavior does not depend strongly 
on p, we well depict only the phase diagrams for K/J = 0, where qualitative 
changes occur when p is varied. 

In Figs. 1 and 2 we show the kT/zJ x A/zJ phase diagrams for K/J = 
and for some values of p. For r = 0.1, the reentrant behavior observed for 
p — 1 is lost and a tricritical point is introduced for p = 3. This latter curve is 
qualitatively equivalent to the phase diagram for r close to zero. For r = 0.3 
a reentrant behavior is obtained for high enough p (curve(c)); a critical end 
point (E), which was not present for p — 1 and p = 3, is also observed. The 
first order transition curve for p = 7 extends to negative values of A/zJ, until 
it reaches a line of continuous transitions at a tricritical point (not depicted 
in the Fig. 2). 

3 Conclusion 

Using a mean-field approximation, we calculate the phase diagrams for the 
random crystal field DBEG model, varying the degree of disorder, r, and 
the degeneracy factor p. For high values of K/ J the phase diagrams are al- 
most insensitive to the value of p, while for K/J = the dependence on p is 
stronger. For the latter value of K/J and r = 0.3 the phase diagram changes 
qualitatively when p is increased enough: while for p = 1, for instance, the 
order-disorder transition is always continuous and the first-order transition 
inside the ordered phase ends in a critical point, for p = 7 a line of first or- 
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der transition separates the ordered and disordered phases at low values of 
A/zJ. A reentrant behavior is obtained for high p. For K — and r = 0.1, 
the degeneracy reintroduces a tricritical point and eliminates the reentrant 
behavior. 

Some of the features found in this work may be artifact of the simple mean- 
field approximation we used and should not hold for a low-dimensional system. 
Preliminary results for the square lattice, using a real-space renormalization- 
group approach, show a different behavior: for any infinitesimal amount of 
disorder the first order transition is supressed, independent of p. 
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FIGURE CAPTIONS 



Figure 1: Phase diagram of the random-field DBEG model, for K/J = and 
r = 0.1. Dashed (continuous) lines indicate first-order (continuous) transition, 
k is the Boltzman constant, T is the temperature and z is the coordination 
numer of the lattice. 01(02) stands for the ordered phase with m — l(m = r) 
at T = 0, D stands for the disordered phase, C stands for critical points, E 
stands for critical end points, and TC stands for tricritical points. Finally, 
(a) stands for p = 1, while (b) stands for p = 3. The diagram for r = is 
qualitatively the same as the curve (b), except for the 02 phase, which is no 
longer present. 

Figure 2: Phase diagram of the random-field DBEG model, for K/J = and 
r = 0.3. Same notation as in previous figure; (c) stands for p = 7. Note the 
qualitative change for p > 3: a critical end point is present for p = 7, as well as 
reentrant behavior, with order-disorder-order transitions as A/zJ is increased 
from small values. 
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